Rolling friction for hard cylinder and sphere on viscoelastic solid 
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We calculate the friction force acting on a hard cylinder or spherical ball rolling on a flat surface 
of a viscoelastic solid. The rolling friction coefficient depends non-linearly on the normal load and 
the rolling velocity. For a cylinder rolling on a viscoelastic solid characterized by a single relaxation 
time Hunter has obtained an exact result for the rolling friction, and our result is in very good 
agreement with his result for this limiting case. The theoretical results are also in good agreement 
with experiments of Greenwood and Tabor. We suggest that measurements of rolling friction over 
a wide range of rolling velocities and temperatures may constitute an useful way to determine the 
viscoelastic modulus of rubber-like materials. 



1 Introduction 

Rubber friction is a topic of huge practical importance, 
e.g., for tires, rubber seals, wiper blades, conveyor belts 
and syringes [l- 1A\. Many experiments have been per- 
formed with a hard spherical ball rolling on a flat rubber 
substrate 1|1 17, 18 1 . Nearly the same friction force is 
observed during sliding as during rolling, assuming that 
the interface is lubricated and that the sliding velocity 
and fluid viscosity are such that a thin lubrication film is 
formed with a thickness much smaller than the indenta- 
tion depth of the ball, but larger than the amplitude of 
the roughness on the surfaces [15(. The results of rolling 
friction experiments have often been analyzed using a 
very simple model of Greenwood and Tabor [l5j|. which 
however contains a (unknown) factor a, which represent 
the fraction of the input elastic energy lost as a result 
of the internal friction of the rubber. In this paper we 
present a very simple theory for the friction force acting 
on a hard cylinder or spherical ball rolling on a flat rub- 
ber surface. For a cylinder rolling on a viscoelastic solid 
characterized by a single relaxation time Hunter [lfjj has 
obtained an exact result for the rolling friction, and our 
result is in very good agreement with his result for this 
limiting case. 

2 Theory 

Using the theory of elasticity (assuming an isotropic 
viscoelastic medium) , one can calculate the displacement 
field Ui on the surface z = in response to the surface 
stress distributions Cj = a^. Let us define the Fourier 
transform 
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d 2 x dt Ui(x, t)e 
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and similar for <Xi(q, u>). Here x = (x, y) and q = (q x , q y ) 
are two-dimensional vectors. In Ref. [3[ we have shown 
that 

Ui(q,cd) = My(q,a;)<Tj(q,a>) 
or, in matrix form, 

u(q,w) = M(q,u>)a(q,u>) 



where the matrix M is given in Ref. 

We now assume that | Vw z (x) |< 1 and that the sur- 
face stress only acts in the z-direction so that 



u z (q,w) = M zz (q,u)a z (q,u). 



(1) 



Since in the present case u> is of order vq (where v is the 
sliding or rolling velocity) we get uj/crq — v/ct « 1 in 
most cases of practical interest, where ct is the transverse 
sound velocity in the rubber. In this case (see Ref. Q): 
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It is interesting to note that if, instead of assuming that 
the surface stress act in the z-direction, we assume that 
the displacement u point along the z-direction, then 

cr 2 (q,cj) = (M -1 )^ (q,u)u z (q,u) 

where in the limit cj/opq << 1, 
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which differ from (2) only with respect to a factor 
4(1 - ^) 2 /(3 - 4u). For rubber-like materials (u « 0.5) 
this factor is of order unity. Hence, practically identical 
results are obtained independently of whether one as- 
sumes that the interfacial stress or displacement vector 
is perpendicular to the nominal contact surface. In re- 
ality, neither of these two assumptions hold strictly, but 
the result above indicate that the theory is not sensitive 
to this approximation. 
Now, assume that 



er z (x,t) = er z (x - vi) 



then 



*,(q,w) = 7^3 J d 2 xdt o-,(x-vt) e-^ x -"') 
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where 



«r«(q) = 7A2 y d2x °"*( x ) e iq ' : 



(27T) 

If Ff denote the friction force then the energy dissipated 
during the time period to equals 



AE = F{vt . 
But this energy can also be written as 

AE = I d 2 xdi u ■ cr 

3 / j2 



(4) 



= (2n) J d qdu (— iu>)u(q, uS) ■ a(— q, — uS) (5) 

where oj = v ■ q. Substituting (1) in (5) and using (3) 
and that 

[<5( w -q-v)] 2 = (to/27r) 5(w-q-v), 

gives 



AE = (2tt)% / d 2 q (-iu) 



xM zz (-q, -u)a z (q)a z (-q). 

Comparing this expression with (4) gives the friction 
force 



(27T) 



■y*d 2 g (-iu) M zz (-q, -u)a z (q)a z (-q). 
Since Ff is real and since a z (— q) = cr*(q) we get 

Ff = I d2q wImM ^(~ q ' - w )Mq)l 2 - ( 6 °) 



Using (2) we can also write 
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d q — Im 



q E cS (u) 
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where F c g — E/(l — v 2 ). In principle, depends on fre- 
quency but the factor 1/(1 — v 2 ) varies from 4/3 « 1.33 
for v = 0.5 (rubbery region) to ~ 1.19 for v = 0.4 (glassy 
region) and we can neglect the weak dependence on fre- 
quency. Within the assumptions given above, equation 
(6) is exact. Note that even if we use a z (q) calculated to 
zero order in tan<5, the friction force (6) will be correct 
to linear order in tanJ. We also note that the present 
approach is very general and flexible. For example, in- 
stead of a semi-infinite solid as assumed above one may 
be interested in a thin viscoelastic film on a hard flat sub- 
strate. In this case, (assuming slip-boundary conditions) 



the M(q, o;)-function which enter in (6a) is determined 
by@,Q 
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4i/)cosh(2gd) + 2(qd) 2 - 4i/(3 - 2v) + 5 



(3 - 4^)sinh(2<7d) - 2qd 



Note that as d — > 00, S — > 1 and the present result 
reduces to (2). Substituting this in (6a) gives the rolling 
friction for a sphere (or cylinder) on a thin rubber film 
adsorbed to a hard flat substrate. We now apply the 
theory to (a) a rigid cylinder and (b) a rigid sphere rolling 
on a semi-infinite viscoelastic solid. 

2.1 Cylinder 

Consider a hard cylinder (radius R and length L y >> 
R) rolling on a viscoelastic solid. The same result is ob- 
tained during sliding if one assume lubricated contact and 
if one can neglect the viscous energy dissipation in the 
lubrication film. As discussed above, when calculating 
the friction force to linear order in tani5 we can neglect 
dissipation when calculating the contact pressure a z (x) 
and assume that the stress is of the Hertz form. Thus if 
we introduce a coordinate system with the j/-axis parallel 
to the cylinder axis and with the origin of the x-axis in 
the middle of the contact area (of width 2a), then the 
contact stress for —a < x < a: 
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where /n = F^/L y the load per unit length of the cylin- 
der. The half-width of the contact area in the Hertz 
contact theory is 
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where we take E e g to be \E c g (u)\ with uj — q • v. Now, 
note that 
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FIG. 1: The friction coefficient (times the radius R of the 
cylinder and divided by the half- width ao of the static contact 
region) as a function of vr/ao, where v is the rolling velocity, 
r the rubber viscoelastic relaxation time (for Ei/Eq = 10). 
We compare the exact result (blue curve) of Hunter QJ] with 
the prediction of Eq. (11) (green curve). 



Using (7) and (9) gives 



<*(<v 
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S(q y )Ji(q x a) 



(10) 



Note that this expression satisfies 

(2^) 2 ( r(q = 0) = J d 2 x a z {x) = F N 

where we have used that S(q y = 0) = L v /2tt. Substitut- 
ing (10) in (6b) dividing by Fr, and using that 

[5(q y )} 2 = (L y /2n)S(q y ), 

gives the friction coefficient: 
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;Jl{q x a) (11) 



Let us now assume the simplest possible viscoelastic 
modulus characterized by a single relaxation time r: 



Exil-iwr) 

Ei/Eq — iLOT 



(12) 



where E\ / Eq is the ratio between the high frequency and 
low frequency modulus. In Fig. Q]we show the friction co- 
efficient (times the radius R of the cylinder and divided 
by the half-width ao = a(v = 0) of the static contact 
region) as a function of vr/ao, where v is the rolling ve- 
locity. We have assumed Ei/Eq = 10. We compare the 
exact result (blue curve) of Hunter [3] (the same result 
was obtained by Goryacheva 22] ) with the prediction of 
Eq. (11) (green curve). Note that some distance away 
from the maximum the agreement between the two curve 



is perfect. This is expected because these regions corre- 
spond to small tan<5 where Eq. (11) should be essen- 
tially exact. Close to the maximum a small difference 
occur between the two curves, but from a practical point 
of view this is not important, since real rubber exhibit 
some non-linearity making any linear viscoelasticity the- 
ory only approximately valid anyhow. 

At this point we empathize that (6) is basically exact, 
and if one could calculate the contact pressure a z (x) [or 
rather the Fourier Transform <J z (q x )] exactly, then (6) 
should give the exact result, e.g., the result of Hunter for 
the cylinder case, and assuming the simple viscoelastic 
modulus (12). Note that the contact pressure will not 
be symmetric around the midpoint when the dissipation 
in the rubber is included in the analysis. Still, to lin- 
ear order in tanJ one can neglect this effect [since (6) is 
explicitly already linear in tan<5], and the analysis above 
shows that this is a remarkable accurate approximation, 
which is an interesting result in its own right, and makes 
it possible to apply the theory to a wide set of problems. 

2.2 Sphere 

Consider a hard spherical (radius R) ball rolling on a 
viscoelastic solid. As discussed above, when calculating 
the friction force to linear order in tani5 we can neglect 
dissipation when calculating the contact pressure er z (x) 
and assume that the stress is of the Hertz form: 
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(13) 



where r — |x| is the distance from the center of the con- 
tact area and F^ the load. The radius of the contact area 
in the Hertz contact theory is 



V 4£ cff J 



(14) 



where we take E e g to be \E c g(uj)\ with uj = q • v. Now, 
note that 
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FIG. 2: The friction coefficient as a function of vr/d, where 
v is the rolling velocity, r the rubber viscoelastic relaxation 
time, and d — ao for the cylinder (upper curve) or d = ro 
for the sphere (lower curve), rolling on a rubber substrate 
described by a simple viscoelastic model (see Eq. (12)). Here 
ao is half the width of the cylinder-substrate contact area and 
ro the radius of the contact region for the sphere, both in the 
limit of vanishing rolling velocity. In the calculation the radius 
of the sphere and of the cylinder are both R = 1 cm and the 
load was chosen so that the average contact pressure for the 
stationary contact was the same (equal to about 44 kPa). For 
E = 1 MPa and Ei/E = 10. 

= 2nrl(qr c y 3 [sin(gr c ) - qr c cos(qr c )] (15) 
Using (13) and (15) gives 

a ^ = ( g7 ! )3 [ sin (9 r c) ~ qr c cos(qr c )} (16) 

Note that this expression satisfies 

(27r) 2 cr(q = 0) = / d 2 x er z (x) = F N 

Substituting (16) in (6b) and dividing by i"N gives the 
friction coefficient: 

9F N [°° f 2 \ 1 
(j, = —r / dqq / d(j> cos</> Ira.-— — 

1 2 

x r s 6 [sin(gr c ) - gr c cos(gr c )] (17) 

Eq. (17) is very general with no restriction on the vis- 
coelastic properties of the rubber or on the sliding veloc- 
ity v (assuming that v is small enough that the effect of 
frictional heating is negligible, and assuming v << op). 

In Fig. [2] wc show the rolling friction coefficient as 
a function of vt/tq for a sphere, and as a function of 
vr/ao for a cylinder. The ball has the radius R = 1 cm 
and is squeezed against a rubber surface with the load 
F/v = 0.15 N giving a static contact area with the radius 
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FIG. 3: The maximum (as a function of the velocity v) of 
the rolling friction coefficients for a sphere and a cylinder, as 
a function of the average (static) contact pressure. For the 
same parameters as in Fig. [2] 

r = r c (v = 0) = 0.10 cm. The average contact pressure 
Pa = Fv/( 7rr o) = 44.3 kPa is the same as for the cylin- 
der case shown in the same figure. The cylinder has the 
radius R = 1 cm and is squeezed against a rubber sur- 
face with the load /jv = 100 N/m giving a static contact 
area with the width ao = 0.11 cm. The average contact 
pressure is p a = /jv/(2ao) = 44.3 kPa. In the calculation 
we have used the viscoelastic modulus E c g(uj) given in 
Eq. (12). Note that at very low sliding velocities where 
the (average) contact pressure are the same the rolling 
frictions are nearly identical. 

In figure |3] we show the maximum (as a function of the 
velocity v) of the rolling friction coefficients for a sphere 
and a cylinder, as a function of the average pressure in 
the static contact area. Note that fj, varies nearly linearly 
with the (static) contact pressure, and using (8) and (14) 

1 /2 

it follows that the rolling friction coefficient /i ~ / N for 
the cylinder and /x ~ F^/ 3 for the sphere. 

Let us now consider the limiting case when the rolling 
velocity v is so small that only the low-frequency vis- 
coelastic modulus is relevant. We also assume that 
that E(uj) is given by (12). If wt << 1 for all rel- 
evant frequencies, we get Re_E(u;) w E(0) — Eq and 
IhlE(u;) w -utE^I-Eq/Ex). We also get \E{uj)\ w Bo- 
Substituting these results in (17) gives after some simpli- 
fications 

^ 2 E V Eh.) r V > 

where p a = Fn/^titq) is the (average) static contact pres- 
sure, and where 

/>oo 

/ = / d£ C 4 (si< - £cos0 2 ~ 0.52 
Jo 

Note that fi ~ p a with is consistent with Fig. [3J If 
we define the frequency lj = v/tq and note that in the 



present case t&n5 = (1 — Eo/E\)ujt we can write 



— — t&nS. 
2 E 



Since 91/2 w 2.34 we get 



2.34-p-tan<5. 



(19) 



In Fig. 2] we compare this asymptotic (low velocity) re- 
sult with the full theory [Eq. (17)]. For the limiting case 
studied above Greenwood and Tabor have shown that 



97T Pa 

64 E 



a, 



(20) 



where a is the fraction of the input elastic energy lost as 
a result of the internal damping of the rubber. Thus we 
can write 



a « 5.3tan5. 



(21) 



Greenwood and Tabor analyzed rolling (and lubricated 
sliding) friction data using (20). They obtained the best 
fit to the experimental data by using an a which was 
almost a factor of two larger than obtained from cyclic 
(low frequency) simple uniaxial loading-unloading mea- 
surements of rubber strips, where the energy loss due 
to hysteresis was about 0.35 of the maximal elastic en- 
ergy of deformation. For a linear viscoelastic solid the 
latter is given by irtand (see Appendix A). Our theory 
predict that a (see (21)) is about 1.7 times larger than 
7rtan<5, and using this result gives very good agreement 
with the measured data presented in Ref. [15j . The rea- 



son for why rolling friction experiments result in an a- 
parameter which is larger than expected from uni-axial 
tension tests is related to the very different nature of the 
time-dependent deformations: a rubber volume element 
below a rolling sphere (or cylinder) undergoes, for some 
fraction of the interaction time, time-dependent defor- 
mations where the total elastic energy is nearly constant 
while the stress directions are chan ging and energy is be- 
ing lost (see discussion in Ref. [l5, 16|). 



3 Rolling friction on rubber 

We have calculated the rolling friction when a hard 
cylinder with radius R = 1 cm is rolling on a Styrene- 
Butadiene (SB) copolymer rubber surface at room tem- 
perature. The viscoelastic modulus of the rubber has 
been measured, and in Fig. [5] we show the loss tangent 
tan<5 = ImE(u})/ReE(uj) at T = 20 C, as a function of 
frequency w for an unfilled and filled SB-rubber. Note 
that tan<5 for the filled SB-rubber has a tail extending 
to very low frequencies. This correspond to relaxation 
processes in the rubber with very long relaxation times 
(see Fig. B.2 in 0]). This effect result from the increase 
in the activation barrier for polymer rearrangement pro- 
cesses for rubber molecules bound to (or close to) filler 
particles (in this case carbon particles). 
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FIG. 4: The friction coefficient as a function of vr/ro for the 
sphere. The full line is using the full theory [Eq. (17)] and 
the dashed line using the asymptotic, low velocity, result (19). 
For the same parameters as in Fig. [2] 
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FIG. 5: The loss tangent tan<5 = lmE(uj)/ReE(uj) at T = 
20 C, as a function of frequency lo for an unfilled and filled 
SB rubber. 



The rolling friction coefficient is shown in Fig. [5] when 
the normal force per unit length is 10 N/cm. The rolling 
friction exhibit a very strong temperature dependence 
given by the WLF shift factor ar'- reducing the temper- 
ature by 10 C result in a shift of the rolling friction curve 
towards lower velocities by approximately one decade. 
At low frequencies (or low rolling velocities) the unfilled 
SB rubber is elastically much softer than the filled rub- 
ber. Thus, the half width ao of the static Hertz contact 
region is 1.4 mm for the unfilled rubber but only 0.6 mm 
for the filled rubber. At high rolling velocities these val- 
ues becomes much smaller owing to the stiffening of the 
effective clastic modulus at high frequencies. Note also 
that the rolling friction for filled SB rubber has a tail 
extending towards lower rolling velocities, which has the 
same origin as the tail in tan<5 towards lower frequencies, 
involving rubber molecules bound to the filler particles. 
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FIG. 6: The rolling friction coefficient for a hard cylinder 
with radius R = 1cm rolling on unfilled and filled SB rubber 
at T = 20 C. The normal force on the cylinder per unit length 
is 10 N/cm. 

4 Discussion 

Consider a rolling cylinder (or ball) in a reference frame 
where the center-of-mass velocity v vanish. The rolling 
friction gives a moment FfR at the center of the cylin- 
der. Since the shear stress at the interface is assumed to 
vanish, this moment must arise from a contact pressure 
a z (x) which is asymmetric around the center line of the 
cylinder. That is, 

FfR — L y J dx xa z (x). 

Thus the calculation of the rolling resistance gives infor- 
mation about how the centroid of the contact pressure 
shift away from the symmetry line x = due to the de- 
layed response of the rubber caused by the internal fric- 
tion of the rubber (characterized by the relaxation time 
r in (12)). This information could in principle be used to 
improve the theoretical treatment given above, but since 
the simple theory we use seems to be very accurate, we 
will not consider this point further here. 

The study presented above assumes that the adhe- 
sional interaction between the rubber and the hard ball 
can be neglected. Most rubber of engineering interest are 
strongly cross-linked and have fillers making them rela- 
tive stiff and reducing the role of adhesion. Adhesion can 
also be removed (or reduced) by lubricating the interface 
or by adsorbing small inert solid particles on the rub- 
ber surface, e.g., talc. In these latter cases, the rolling 
and sliding friction may be nearly the same as indeed 
observed in some experiments [lij] . For very soft rubber 
and for clean surfaces, during rolling or sliding an open- 
ing (and a closing) crack will propagate at the interface, 
and associated with this may be very strong energy dissi- 
pation, which may dominate the rolling or sliding friction 
(see, e.g., Ref. [21|.). We note finally that measurements 
of the rolling friction as a function of velocity v and tem- 



perature T may be a very useful way of determining the 
viscoelastic modulus E{uS), see Appendix B. 

5 Summary and conclusion 

We have presented a very general and flexible approach 
to calculate the rolling resistance of hard objects on vis- 
coelastic solids. The theory can be applied to materi- 
als with arbitrary (e.g., measured) viscoelastic modulus 
E(uj). The theory can be applied to both spheres and 
cylinders rolling on semi-infinite viscoelastic solids, or on 
a thin viscoelastic film adsorbed on a rigid flat substrate, 
or even more complex situations for which the M(q, oj)- 
function can be calculated. For a cylinder rolling on a 
viscoelastic solid characterized by a single relaxation time 
r, Hunter has obtained an exact solution for the rolling 
resistance. We have shown that for this limiting case our 
theory gives almost the same as the result as obtained 
by Hunter. We have compared the rolling resistance of 
a sphere with that of a cylinder under similar circum- 
stances. Measurements of the rolling friction as a func- 
tion of velocity and temperature may be a very useful 
way of determining the viscoelastic modulus. 
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Appendix A 

Consider a strip of rubber in (slow) loading-unloading 
and assume that the strain e = eosin(wi), < t < to with 
uito = ir. The energy dissipation (per unit volume) 

t^diss = / dt ere 

JO 

= [° dt f£ (E(Lu)e 1 ^ - E(-u)e- iut ) 
Jo 2i 

x^w (e lut + e~ iut ) = legwt Im£(w) 

We are considering very low frequencies so that the max- 
imum elastic energy 



Thus we get 



t/di SS lmE(u) . 
— — = 7r „ „, N ps 7rtano 



ReE(u) 
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Appendix B 

The viscoelastic modulus E(u>) is a complex quantity 
and at first one may think that it impossible to de- 
termine both ReE(ui) and IuiE(uj) from a knowledge 
of a single function n(v) which depends on E(uj) = 
KeE(uj) + ilmE(uj). However, E(w) is (assumed to be) 
a causal linear response function so that JmE(ui) can in 
fact be obtained from Re.E'(w) using a Kramer-Kronig 
relation. Thus we have one known function fJ,(v) and 
only one unknown function, e.g, KeE(uj), and E(uj) can 
in principle be determined uniquely from fJ,(v). The best 
way of doing this is to note that E(w) can be written as 



POO 

E(uj) = E x - / dr- 



H(r) 



{Bl) 



In numerical calculations one may discretize the relation 
(Bl) and it is usually enough to include n « 15 relaxation 
times Tfc with Tk+i ~ lOrfe. Thus we expand E(uj) on the 
form 



k=l 



lOJT k 



(B2) 



where the Hk and the relaxation times are real positive 
quantities. Next, form the quantity 

V = y^/Aneasfci) ~ Mtheory^i; H X , H n )} 2 gi (S3) 



where <?, are suitable chosen weight coefficients (e.g., 
gi = 1) and where Vi are the velocities for which the 
rolling friction has been measured. In (B3) we have in- 
dicated that the theory expression for the rolling friction 
coefficient depends on the numbers Hk- The determina- 
tion of Hk is a problem in multidimensional minimization 
of V, and can be performed using different methods, e.g., 
the Monte Carlo method or the Amoeba method, see Ref. 

M. 
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